Infinite-randomness quantum critical points induced by dissipation 
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We develop a strong-disorder renormalization group to study quantum phase transitions with 
continuous 0{N) symmetry order parameters under the influence of both quenched disorder and 
dissipation. For Ohmic dissipation, as reahzed in Hertz' theory of the itinerant antiferromagnetic 
transition or in the superconductor-metal transition in nanowires, we flnd the transition to be 
governed by an exotic infinite-randomness fixed point in the same universality class as the (dissipa- 
tionless) random transverse-field Ising model. We determine the critical behavior and calculate key 
observables at the transition and in the associated quantum Griffiths phase. We also briefly discuss 
the cases of superohmic and subohmic dissipation. 
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I. INTRODUCTION 

In recent years, it has become clear that quenched 
disorder, i.e., impurities, defects or other types of im- 
perfections can significantly modify the low-temperature 
behavior of quantum many-particle systems. At zero- 
temperature quantum phase transitions, the interplay 
between large-scale quantum fluctuations and random 
fluctuations leads to much more dramatic disorder ef- 
fects than at classical thermal phase transitions, resulting 
in various exotic phenomena such as quantum Griffiths 
effects^'^''^, non-power-law dynamical scaling^'^, or even 
smeared phase transitions.®''' A recent review of some of 
these phenomena can be found in Ref. 8. 

The quantum phase transitions in random transverse- 
field Ising magnets are among the most striking examples 
of such behavior. Utilizing a real-space renormalization 
group (RG) technique due to Ma, Dasgupta and Hu,^'^" 
now known as the strong-disorder RG (see Ref. 11 for a 
review), Fisher^'^ showed that the one-dimensional ran- 
dom transverse-field Ising chain features an unconven- 
tional infinite-randomness critical point with ultraslow 
activated rather than power-law dynamical scaling. It 
is accompanied by strong power-law quantum Griffiths 
effects in the vicinity of the transition. While it was ini- 
tially suspected that this scenario is special to one space 
dimension, Motrunich et al}'^ showed that the random 
transverse-field Ising models in two and three dimensions 
also display infinite-randomness critical points. 

A dissipative environment further hampers the dynam- 
ics. In the experimentally relevant case of Ohmic damp- 
ing, the large locally ordered droplets that are normally 
responsible for quantum Griffiths effects completely cease 
to tunnel. ^'^'^"''^^ Instead, they develop static order in- 
dependently from the bulk system which destroys the 
sharp quantum phase transition by smearing.®'^ A sim- 
ilar freezing of locally ordered droplets also occurs close 
to a quantum percolation transition. -'^^'-'^^ 

The above behavior of Ising order parameters must be 
contrasted with that of continuous 0{N) symmetry or- 
der parameters. While the bulk ground state phases of 



one-dimensional Heisenberg random quantum spin chains 
are governed by infinite-randomness fixed points, 
higher-dimensional random quantum Heisenberg sys- 
tems have more conventional ground states, ^^'^^ and 
their quantum phase transitions are governed by con- 
ventional critical points. ^^'^^'^'^'^^ As in the Ising case, 
adding Ohmic dissipation hampers the dynamics of 0{N) 
symmetric order parameters. Vojta and Schmalian^^ 
showed that the "energy gap" of large locally ordered 
droplets is exponentially small in their volume lead- 
ing to power-law quantum Griffiths effects analogous to 
those in the dissipationless random transverse-field Ising 
model. This analogy suggests the important question of 
whether Ohmic dissipation can induce an unconventional 
infinite-randomness quantum critical point for a contin- 
uous 0{N) symmetry order parameter. 

In addition to its significance for the theory of quantum 
phase transitions, this question also has important ex- 
perimental applications. Consider the antiferromagnetic 
quantum phase transition of itinerant electrons (as ob- 
served, e.g., in heavy fermion compounds^®'^^). Within 
the standard Hertz-Millis spin-fluctuation theory^®'^^, it 
is described by an 0(3) Landau-Ginzburg- Wilson (LGW) 
order parameter field theory with Ohmic dissipation. 
The properties of this transition have been a long- 
standing unsolved problem; and given the fact that most 
experimental systems are rather dirty, studying the ef- 
fects of disorder on the Hertz-Millis theory is of prime 
interest. 

A second potential application is provided by the 
pair breaking superconductor-metal quantum phase 
transitions in nanowires. '^'^ It can be described by 
a one-dimensional 0{2) LGW theory with Ohmic 
dissipation. ■^^'■^^'•^'^ There is experimental evidence that 
the pair breaking in this systems is caused by sur- 
face magnetic impurities which necessarily also introduce 
quenched disorder. 

In this paper, we investigate the quantum phase tran- 
sition of a continuous-symmetry 0{N) order parameter 
under the combined influences of both quenched disor- 
der and Ohmic dissipation. To this end we develop a 
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strong-disorder RG suitable for this problem. The paper 
is organized as follows: In Sec. II we introduce our model, 
a dissipative LGW order parameter field theory. In Sec. 
Ill we implement the strong-disorder RG in the large- 
limit and relate it to that of the random transverse-field 
Ising model. We also summarize the solution and the 
resulting critical behavior. In Sec. IV, we calculate key 
observables close to the transition while Sec. V deals with 
the case of non-Ohmic damping. We also show that our 
results do not rely on the large- TV limit. We conclude in 
Sec. VI. 

A short account of part of this work has already been 
published in Ref. 34. 

II. ORDER PARAMETER FIELD THEORY 

We start from a quantum Landau-Ginzburg- Wilson 
(LGW) free energy functional for an A^-component {N > 
1) vector order parameter tp = . . . ,(Pn) in d space 
dimensions. For the above-mentioned itinerant antifer- 
romagnetic quantum phase transition, generically d = 3 
and = 3 while for the superconductor-metal tran- 
sition in nanowires, d = 1 and N = 2. The LGW 
free energy can be derived from an appropriate micro- 
scopic Hamiltonian of disordered electrons using stan- 
dard methods^*'^^'^^ (for a critical discussion of this ap- 
proach see Ref. 36). In the absence of quenched disorder, 
the action of our LGW theory reads 

S = J dydx ip{x)T{x,y)ip{y) + ^ J 'f'^(^) ' (1) 

where x = (x, t) comprises imaginary time r and position 

X, J dx = J dx Jq^^ dr, and u is the standard quartic 
coefRcient. r(a;, y) denotes the bare inverse propagator 
(two-point vertex) whose Fourier transform reads 

r(q,a;„) =r + eo'q'+7o|c^„|'^'" • (2) 

Here, r is the bare distance from criticality (the bare 
gap), ^0 is a microscopic length scale, and ll}„ is a Matsub- 
ara frequency. The nonanalytic frequency dependence of 
r(q, w„) is caused by the coupling of the order param- 
eter to a dissipative bath. We are mostly interested in 
the case of overdamped (Ohmic) dynamics correspond- 
ing to a value of zq = 2. However, to demonstrate the 
special role of zq = 2, we also consider different values of 
Zq. The damping coefRcient 70 depends on the coupling 
of the order parameter to the dissipative bath and the 
spectral density of the bath modes. 

In the presence of quenched disorder, the functional 
form of the order-parameter field theory (1) does not 
change qualitatively, but the distance from criticality r 
becomes a random function of spatial position. Analo- 
gously, disorder appears in ^o,7o and u. 

Let us briefly comment on possible modifications of the 
two-point vertex (2) by mode-coupling effects. For the 
itinerant ferromagnetic quantum phase transition"^^'"^*''^^ 



and the superconductor-metal transition without mag- 
netic impurities, ''^"'■^^ the coupling between the order pa- 
rameter fluctuations and the soft particle-hole excitations 
of the metal leads to a long-range interaction in space, 
represented by a nonanalytic g-dependence instead of the 
simple term. In contrast, in our examples, the term 
remains the leading term because the relevant modes are 
either gapped (for the superconductor-metal transition 
due to magnetic impurities) or couple too weakly to the 
order parameter (in the case of the itinerant antiferro- 
magnetic transition). 

Our goal is the application of the real-space-based 
strong-disorder RG. We therefore need to discretize the 
continuum action (1) in space (but not in imaginary time) 
by defining discrete coordinates Xj and rotor variables 
LPjiyT). These rotors are coarse-grained rather than mi- 
croscopic variables, they represent the average order pa- 
rameter in a volume W large compared to the micro- 
scopic scale ^0 but small compared to the true correlation 
length ^, i.e., (^^(t) j^y dyip{xj + y, r). 

For simplicity, we first consider the large-iV limit of 
our LGW theory. This will allow us to perform all of 
the following calculations explicitly. We will later show 
that the RG fixed point is the same for all iV > 1. The 
resulting discrete large- action reads 

S = T^^(r, + A,+7.K|'/'°) |</>.(a;„)|' 

-T^^<f>i{-UJn)Jij<l>j{i^n) , (3) 

where rt, ji > and the nearest-neighbor interac- 
tions Jij > are random quantities, and (f)j{iUn) = 
</5j(T)e*'^"'^dT is the Fourier transform of the rotor 
variable. The Lagrange multipliers enforce the large- 
A^ constraints = ^ ^'^^ each order parameter 

component ip"^^ at each site «; they have to be deter- 
mined self-consistently. The renormalized local distance 
from criticality at site i is given by = + A^. In the 
disordered phase, all > 0. 



III. STRONG-DISORDER RENORMALIZATION 
GROUP 

The idea of the strong-disorder (Ma-Dasgupta-Hu) 
j^q4,5,9,io consists in the successive decimation of local 
high-energy degrees of freedom. It relies on the disorder 
distributions being broad and becomes exact in the limit 
of infinitely broad distributions. For now, we assume 
that our distributions are sufficiently broad, we will later 
show that their widths diverge at the critical RG fixed 
point, justifying the method. 
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A. Single-cluster solution 



1. Decimating a site 



Let us start by considering a single rotor variable 
(i.e.; a single cluster) with the action 

S-ci = T^(r + A + 7|a.„pA-") |,^K)I' ■ (4) 

The value of the Lagrange multiplier A is determined by 
the length constraint 
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r + A + 7|u;„|2/^o 



(5) 



At zero temperature, the Matsubara sum can be turned 
into an integral, resulting in 



^ ~ 7^ ( du) ; — , 



(6) 



where e = r + A is the renormalized distance from criti- 
cality. 

To proceed, we now need to distinguish super-Ohmic, 
Ohmic and sub-Ohmic dissipation. In the super-Ohmic 
case, zq < 2, the integral can be carried out straight- 
forwardly giving ((/)^) = ce'^^°~^-'''^7~^''/^ with c being a 
constant. Solving for e yields the relation between the 
gap and the damping constant (i.e., the cluster size) 



(1/7) 



zo/{2-zo) 



(7) 



In contrast, we need to introduce a high-frequency cut- 
off A to carry out the constraint integral in the Ohmic 
case Zq = 2, giving (0^) = ln(l + jA/e)/n'-f. The result- 
ing dependence of the gap on the damping constant is 
exponential, 



7A/ (e'^'^ - 1) w 7Ae- 



(8) 



signifying that a single Ohmic cluster is marginal, i.e., 
right at the lower critical "dimension" of the problem. 
In the sub-Ohmic case, zo > 2, the single-cluster physics 
changes dramatically: The constraint integral in (6) con- 
verges in the limit e — > 0. Thus, once 7 > 7c = 
^(zo-2)/zo^g/(7r(zo-2)), Eq. (6) does not have a solution 
e > anymore, implying that the rotor has undergone 
a localization phase transition caused by the sub-Ohmic 
dissipation. 



B. Recursion relations 

In our large- iV action (3), the competing independent 
local energies are the gaps and the interactions Jij 
(the damping coefficient 7; and the gap are not in- 
dependent, they are coupled via the large- constraint 
at site i). In the bare theory, the and are inde- 
pendent random variables with probability distributions 
P{J) and i?(e), respectively. Each step of the strong 
disorder RG eliminates one rotor variable by first identi- 
fying the largest local energy Q ~ max(ei, Jij) and then 
decimating the associated high-energy degree of freedom. 



Specifically, if the largest local energy is a gap, say 
62, the corresponding rotor (f)2 is far away from criticality 
and does not contribute to the macroscopic order param- 
eter. However, integrating out its fluctuations generates 
effective interactions between all pairs of sites that cou- 
ple to 02- If the disorder distributions are broad, 62 is 
much larger than all local energies associated with the 
neighboring sites. Thus, 4)2 can be integrated out in per- 
turbation theory with the unperturbed part of the action 
being 



^0 = r^(e2 + 72 k.|'^'")|</'2(a;„)|' 



(9) 



while Si = S — So is the perturbation. Up to 2nd or- 
der in perturbation theory, we only need to consider the 
interaction of 02 with the neighboring sites j, thus 

Si = -T ^ ^2i02(-W„)0j(w„) . (10) 

The partition function can now be written as 



Z 



= Zn 



D[<f>2]l[D[<f>,] 

j¥2 



(11) 



where J D[(j)j] comprises integration over all frequency 
components of , Zq is the partition function associated 
with the action Sq, and (•)o denotes the average w.r.t. Sq. 
The renormalized action S can be calculated in cumulant 
expansion 

S=-He-^%^{Si)o-l[{Sf)o-{Si)l]±... . (12) 

Evaluating the averages, we obtain (5i)o ~ due to 
symmetry and 



E 



£2 +72|w„p/^0 



+2 



Ji2J2]<P*i{^n)4'j{^n) 

£2 +72|t^«p/^° 



(13) 



The first term in the square brackets just gives sublead- 
ing renormalizations of the gaps tj of the neighboring 
sites and can thus be dropped. The second term pro- 
vides the renormalized interactions Jy between all sites 
that used to couple to 02. Their leading low frequency 
behavior is Jy = Ji2J2j/^2 independent of the exponent 
Zo- This term has to be added to the interaction Jy al- 
ready coupling sites i and j, if any. Consequently, the 
final recursion relation for the RG step reads 



Ji^ — Ji' 



Ji2J'. 



i2'J2j 
£2 



(14) 



4 



At the end of the RG step, 02 is dropped from the action. 
Note that the muhiphcative structure of the effective in- 
teraction in (14) is a direct consequence of second order 
perturbation theory. It does not depend on details of the 
model, in particular, it is valid for any zq. 



2. Decimating a bond 

Let us now consider the RG step in the case of the 
largest local energy being an interaction, say J23 cou- 
pling sites 2 and 3. For broad disorder distributions, 
J2'i ^ €2, £3- Thus, the two rotors 4>2 and ^3 are essen- 
tially parallel and can be replaced by a single rotor 02 
which represents the entire cluster comprising 02 and 03 . 
The moment /i2 of the effective rotor, i.e., the number of 
original sites in the cluster is the sum of the moments ^2 
and ^3 of the original rotors, 



M2 = M2 + M3 ■ 



(15) 



To find the renormalizcd gap £2 of the effective rotor, we 
solve exactly the two-site problem involving 02 and 03 
while treating the couplings to all other sites as pertur- 
bations. The two-site action is given by 

5o = ^(r, + A,+7.Kl'^"')l0.K)P 

u)„ i=2,3 

-T^ J23 02(-Wn)'^3(^n) • (16) 

It is subject to the large- A'' length constraints 

da 



1 = 



1 = 



'2)-r 
1) = tY. 



d2d3 - J74 ' 
d2 

d2d3 - J74 ' 



(17) 



with dj = rj + Xj + 7j|a'„p/^". They determine the La- 
grange multipliers A^. (It is important to note that the 
value of ri + Xi in the two-site cluster is different from 
the single-site e,.) 

To integrate out the high-energy mode, we diagonalize 
the quadratic form in (16) separately for each Matsubara 
frequency. The two eigenvalues read 



Ka,b 



d2 + d3± J {d2 



dzY + Jh 



{d2 +d3±J: 



23) 



J23 J23 



(18) 



The corresponding eigenmodes are given by ■0a = a02 
/303 and 'ipb = ~/302 + ^03 with 



d3-d2 + Vid2 - d^r + .P 



{ds 



J 



{ds 



(19) 



2 - 4)2 + J2)2 + J2 



The higher eigenvalue Kb is at least J23 above the 
lower eigenvalue Ka] we thus integrate out the cor- 
responding mode leaving us with the effective action 
S = Ao|'0a(w„)P and a length constraint {xj}"^) = 

{{a<j>2+(ih?)^l. 

We define the renormalizcd rotor variable by rescaling 
02 = i^a/ because wc wish it to fulfill the same 
length constraint (01) = 1 as all other rotor variables. 
Inserting this definition into the (diagonalized) two-site 
action (16) allows us to identify the renormalizcd gap, 
damping constant, and interactions with the neighbors. 



£2 

72 

J2J 



(V'a)(r2+A2+r3-f A3- J) 



= I^{ra){l2+13) , 

= (V^^)'/'(a^2j+/?4j)L„->0 



(20) 



To proceed further, we need explicit results for the La- 
grange multipliers A2 and A3 as well as ("00.) ■ This re- 
quires the solution of the two coupled integral equations 
(17). In the case of Ohmic dissipation, the integrals are 
rational and can be done exactly. In the limit of strong 
disorder, we obtain r2 -|- A2 -|- r'3 -I- A3 — J = 26263/ J and 
aU„^o = /3U„^o = V2/2. Moreover, (0;^) is bounded 
between 1 and 2 and approaches 2 in the asymptotic limit 
^ 0. This leads to the recursion relations 



J: 



2j 



■hj + J: 



3j 



12 = 2 



£263 



J'. 



23 



(21) 



(22) 



implying an additive relation for the renormalizcd damp- 
ing constant. 



72 = 72 + 73 



(23) 



Wc emphasize that the multiplicative form of (22) is not 
independent of the functional form of the action (3). In 
contrast to the recursion relation (14) for the interac- 
tions, the recursion relation (22) for the gaps is special 
to the case of Ohmic dissipation. It is related to the fact 
that the gap e of single cluster depends exponentially 
on the damping constant (and thus on the cluster size), 
e = ^he~'^'^ , as derived in (8). We will come back to this 
point in Sec. V where we discuss the case of non-Ohmic 
dissipation. 

Although the prefactor in Eq. (22) is larger than 1, 
this does not mean that the renormalizcd gap can be- 
come larger than the decimated ones in the weak disor- 
der limit. Using the methods of Ref. 43 we showed that 
the exact value (20) of I2 [calculated within the two-site 
action (16)] is always less than the decimated gaps (£2 
and £3) for all £2,63 < ^2- Therefore, the system flows 
towards the infinite-randomness fixed point for all bare 
disorder strengths, ensuring the internal consistency of 
the RG. 
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The net result of a single RG step is the elimination of 
one rotor and the reduction of the maximum local energy 
n as well as renormalizations of the remaining energies 
and reconnections of the lattice. 

The RG recursion relations (14), (15), (21) and (22) 
completely define the RG procedure. They are identical 
to the corresponding relations for the dissipationless ran- 
dom transverse-field Ising model. "^'^'^^ We thus conclude 
that our system belongs to the same universality class. 
Note, however, that there are some subtle differences in 
the behavior of some observables due to the continuous 
symmetry of the order parameter and the Ohmic damp- 
ing, as will be discussed in Section IV. 

C. RG flow equations and fixed points 

In this subsection, we briefly summarize Fisher's 
solution^^^ of the strong-disorder RG defined by the re- 
cursions (14), (15), (21) and (22) to the extent necessary 
for the purposes of this paper. 

In one space dimension, the RG step docs not change 
the lattice topology because the interactions remain be- 
tween nearest-neighbor sites only, and the e and J remain 
statistically independent. Therefore, the theory can be 
formulated in terms of the probability distributions P{J) 
and i?(e). Fisher derived RG fiow equations for these dis- 
tributions and solved them analytically. There are three 
types of nontrivial fixed points corresponding to the or- 
dered and disordered quantum Griffiths phases and the 
quantum critical point that separates them. 

The most remarkable feature of the critical fixed point 
is that the probability distributions P and R broaden 
without limit under renormalization, even on a logarith- 
mic scale. Using logarithmic variables F = \n{^lj/fl) 
[where flj is of the order of the initial (bare) value of il] , 
( — \n{il/J) and {3 = \n{^l/e) the probability distribu- 
tions 7^(C) and Tl{(3) at the critical fixed point read 

^(C) = ^e-</r, 7^(/3) = ie-z^/i^ . (24) 

The diverging widths of the probability distributions give 
the critical point its name, viz. infinite-randomness crit- 
ical point. They also a posteriori justify the method, be- 
cause the perturbative recursion relations (14)-(22) be- 
come exact in the limit of infinitely broad distributions 
(i.e., approaching the critical point). 

The complete critical behavior can be found by includ- 
ing the moments and lengths of the clusters in the RG 
procedure. It is characterized by three exponents v = 2, 
■0 = 1/2, and = (1 + a/5)/2. The correlation length ex- 
ponent v determines how the average correlation length 
f diverges if one approaches the critical point via 

e~M-^. (25) 

Here r denotes the fully renormalized dimensionless dis- 
tance from criticality which is given by r [ln(e)] — 



[ln(J)] in terms of the bare variables ([•] denotes the dis- 
order average).'*"' 

The tunneling exponent ip controls the dynamical scal- 
ing, i.e., the relation between length scale L and energy 
scale ri. It is of activated rather than power-law type 

ln(r2//rj) - L''' , (26) 

which is a direct consequence of the multiplicative struc- 
ture of the recursions (14) and (22). ip also controls the 
density n^^ of clusters surviving at an energy scale in 
the RG procedure. Its scaling form is given by 

nn{r) = [ln(fi,/0)]-''/'^X„ [r^^^' ln(r!j/0)] , (27) 

with the scaling function behaving as X„(0) = const 
and Xn{y oo) ~ yd/i'g-cdy ^]^Qj.e c is a con- 
stant. As a result, the cluster density decreases as 
~ [ln{D,i at criticality while it behaves as 
nn ^dvQd/z ^Yy^ disordered quantum Griffiths phase 
(r > 0). The dynamical exponent z varies with z ~ r~'^'^ 
in the Griffiths phase. 

Finally, the exponent (p describes how the typical mo- 
ment /in of a surviving cluster depends on the energy 
scale f2. The scaling form of /.j^ reads 

Mo(r) = [Hniin)]^X^ [r-^' ln(r!,/f])] . (28) 

The scaling function behaves as ^^(0) = const and 
X^i{y oo) ~ y^"'^. Thus, at criticality the typical 
moment increases as ~ [In(r2//r2)]'^ while it behaves 
as /io ~ r'"''(*~'^Mn(r2//f2) in the disordered quantum 
Griffiths phase. 

The strong-disorder RG steps discussed in subsection 
IIIB generate effective interactions between sites that 
were previously uncoupled. In dimensions d > 1, this 
changes the lattice connectivity, and it introduces statis- 
tical correlations between the J and e. Therefore, the 
theory cannot be formulated in terms of individual prob- 
ability distributions of these variables, and a closed-form 
analytical solution appears to be impossible. However, 
Motrunich et al}^ numerically implemented the recur- 
sion relations (14)-(22) in two dimensions, keeping track 
of all reconnections of the lattice under the RG. They 
found an infinite randomness critical point very simi- 
lar to that in one dimension. In fact, the critical be- 
havior described in Eqs. (25)-(28) is also valid in two 
dimensions, but with different exponent values. Var- 
ious numerical implementations*^'^^'"'^''*^ of the strong 
disorder RG yielded ip = 0.42... 0.6, = 1.7... 2.5 
and V = 1.07... 1.25. In three dimensions, the RG 
flow towards an inflnite-randomness fixed point has been 
confirmed"'^ , but reliable estimates of the exponent values 
are still missing. 

The strong-disorder RG allows one to identify the 
infinite-randomness fixed point and confirm its stability 
but, strictly, it cannot answer the question of whether or 
not a weakly or moderately disordered system will fiow 
towards this fixed point because, if the disorder is weak. 
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the strong disorder RG step is not very accurate. [An in- 
ternal consistency check^'^ of the RG in the weak disorder 
hmit can be achieved by computing exactly rather than 
perturbatively the renormalized couplings (gaps and in- 
teractions) within the relevant two-site or three-site clus- 
ters, see Sec. IIIB.] For our system, additional insight 
can be gained from the results of a conventional pcrtur- 
bative (replica based) rcnormalization group. Building 
on earlier work,'''* Kirkpatrick and Belitz^^ showed that 
the perturbative RG always takes the system to large 
disorder strength even if the bare disorder is very small. 
Moreover, by taking rare region effects into account in 
an approximate way, Narayanan et aZ.^^'^" showed that 
there is no stable weak-disorder fixed point; instead the 
perturbative RG shows runaway flow towards large disor- 
der. This strongly suggests that our infinite-randomness 
critical point is universal and governs the quantum phase 
transition for all nonzero disorder strength. 



r + \{T) as function of temperature. To this end we solve 
the finite-temperature constraint equation (5) yielding 



6(T) 



e(0) + aT 
T 



(31) 



where e(0) is the zero-temperature value determined by 
the constraint integral (6). In the super-Ohmic and 
Ohmic cases, e(0) is given by Eqs. (7) and (8), respec- 
tively. The constant a is given by a = 7r7 in the Ohmic 
case and a = 2/(2 — zq) in the super-Ohmic case. 

If the rotor variable <f> represents a cluster of moment 
(number of sites) /x, its contribution to the uniform sus- 
ceptibility is proportional to //^ while the contribution to 
the local susceptibility is proportional to fi. By combin- 
ing this with (30) and (31), we obtain the uniform static 
order parameter susceptibility as a function of tempera- 
ture of a cluster of moment /i and distance e from criti- 
cality. 



IV. OBSERVABLES 

The strategy for calculating, within the strong-disorder 
RG, thermodynamic observables such as the susceptibil- 
ity as a function of temperature consists in running the 
RG from the initial energy scale fi/ down ioQ. = T. The 
high-energy degrees of freedom eliminated in this proce- 
dure generally do not make significant contributions to 
the low-energy behavior of observables. At best, they 
change non- universal constants. All clusters surviving at 
energy scale D, = T can be considered to be independent 
because they are coupled by interactions J much smaller 
than T . The desired observable is thus simply the sum of 
independent contributions from the individual surviving 
clusters. Frequency-dependent observables can be deter- 
mined analogously. 



A. Single-cluster results 



Xci(r) 



I^VT [e « T) 



(32) 



The corresponding results for the average local suscepti- 
bility read 



Xcr(r) 



l^/T (e«r) 



(33) 



To calculate the specific heat we also need the total 
energy contribution of a single cluster which behaves as 



ASei-T (e«T). 



(34) 



This is an important difference from the random 
transverse- field Ising model case^ and stems from the fact 
that our rotor variables have an unbounded spectrum. 

We now turn to the dynamical order parameter suscep- 
tibility at zero temperature (focusing on Ohmic dissipa- 
tion). From (30), we obtain in imaginary time formalism 



In order to proceed, we therefore need to calculate the 
relevant observables for single clusters. To do so, we add 
a source term to the single-cluster action (4). It reads 



(29) 



with H{uin) = Jq^^ iJ(T)e*'^"'^dr being the Fourier trans- 
form of the source field conjugate to the order parameter. 
Because the theory defined by Sc\ + Sh is still Gaussian, 
the partition function Zh in the presence of the field can 
be easily evaluated. The dynamic (Matsubara) suscepti- 
bility is then given by 



x(«w„) = 



1 



1 O^IuZh 



(30) 

For the temperature-dependent static susceptibility, we 
set = and find the distance from criticality, e(T) = 



f- + 7|w„| 



(35) 



After Wick rotation iujn ^ w-l-iO to real frequencies, this 
leads to Xc\{'^ + *0) = /^^/(e — i^ijj) implying 



Imxci(w + «0) 



£2 -|- 72,^2 



(36) 



Analogously, the dynamical local susceptibility reads 



-I- 7-^0; 



2,, ,2 



(37) 



B. Summing over all clusters 

We now combine the single-cluster observables summa- 
rized in the last subsection with the strong-disorder RG 
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results for density and moment of the surviving clusters 
given in (27) and (28). We focus on the critical point and 
the disordered Griffiths phase. On the ordered side of the 
transitions, the scenario is dimensionality-dependent be- 
cause in d > 1, an infinite percolating RG cluster forms 
already at a finite energy scale. 

To obtain the uniform static order parameter suscep- 
tibility x(r, T) and the corresponding local susceptibility 
^loc j-j,^ ji-j ^ j.^^ ^]^g pj^Q ^Yie energy scale il = T and 

sum over all surviving clusters. Using (27), (28) and (32), 
we obtain the scaling form 

= ^[HQi/T)]"^-''/^e^[r-'^Hni/T)] (.38) 

with the scaling function given by Q^iy) ~ 

Xn{y)X'^{y). At criticality, r = 0, this leads to x ^ 

[In(r2//T)]^'^"''/'^ /T. In the Griffiths phase we need to 
use the large- argument limit of the scaling function giving 
^ ^ y<i/^-i^di.-H2..V'(i-0)in2(f]^/T). Thus, X shows the 
nonuniversal power-law temperature dependence charac- 
teristic of a quantum Griffiths phase. For z > d, the 
susceptibility actually diverges with T ^ 0. Along the 
same lines, the scaling form of the local susceptibility is 
found to be 

x'°^(r,T) = 1 [In(f7,/T)]'^-'^/'^' e;^- [r-nn(f7,/r)] , 

(39) 

with e^°'=(y) = Xn{y)Xf,{y). This reduces to 

x'°'= - [ln{fli /T)f~'^^'^ /T at criticality and to x'°'' ~ 
j,d/z^i^d^+^i,{i-4,) in(n7/T) in the disordered Griffiths 
phase. 

The scaling form (38) of the susceptibility can also be 
used to infer the shape of the phase boundary close to 
the quantum phase transition. The finite-temperature 
transition corresponds to a singularity in Q^iv) some 
nonzero argument yc- This yields the unusual form ^ 
exp(— const |r| "'"'') shown in Fig. 1. The crossover line 
between the quantum critical and quantum paramagnetic 
regions displays analogous behavior. 

The specific heat C can be found by first adding the 
total energy contributions of all surviving clusters, 

AE{r,T) = TnT{r) 

= T[lnini/T)]-''/^eE[r''^H^i/T)] (40) 

with Qsiy) ~ Xn{y). After taking the temperature 
derivative this gives C - [In(r2//T)] at criticafity 
and C ~ r'^'^T'^/^ in the Griffiths phase. 

To calculate the dependence of the low-temperature 
order parameter m on an external (conjugate) field H, 
we run the RG to the energy scale fijy = ^qH T. All 
decimated clusters have e 2> fiH and do not contribute 
significantly to the order parameter. All surviving clus- 
ters have e ^ iiH and are fully polarized. Summing over 




IRFP 

FIG. 1: (Color online) Temperature-coupling phase diagram 
for Ohmic dissipation. IRFP denotes the infinite-randomness 
critical point. The ordered phase is divided into a conven- 
tional (CO) region and a quantum Griffiths (GO) region. On 
the disordered side of the transition, there is a quantum Grif- 
fiths paramagnet (GPM) followed by a conventional quan- 
tum paramagnet (CPM). The phase boundary (solid) and the 
crossover line (dashed) between the quantum critical (QC) 
region and the quantum paramagnetic regions take unusual 
exponential forms leading to a wide quantum critical region. 
In the classical critical region (CC) close to the phase bound- 
ary classical thermal fluctuations dominate. At sufBciently 
high temperatures (above the dotted dome), the behavior is 
nonuniversal. 

all surviving clusters therefore gives 
m{r,H) = riji„(r)/xo^(r) 

where Qmiv) — Xn{y)X^{y). After resolving the im- 
plicit field dependence caused by the moment in the 
definition of the energy scale fin = fJ-nH, we find 
m - [ln(0//iJ)]'^-''/'^ (with double-logarithmic cor- 
rections) at criticality, r = 0. For r > 0, we ob- 
tain a nonuniversal power-law field dependence, m ~ 
jjd/z^d^+^i,ii-4,){i+d/z) [in{n/H)f+''/\ characteristic of 
a quantum Griffiths phase. 

Finally, to find the zero-temperature dynamic suscep- 
tibility X external frequency lu, we run the RG to 
the energy scale ftui = 70; = ^q^^lo (7 ^ /i7o is the 
effective damping constant of a cluster of moment /i). 
All decimated clusters (having e ^ 70;) only make neg- 
ligible contributions to x- The surviving clusters have 
e ^ 70; which simplifies (36) to Imxci(^) = l^^lla^- In 
the same limit, the local dynamic susceptibility reads 
lmx^°{'{uj) = I/70W. Using (27) and (28) we now sum 
over all surviving clusters to obtain the scaling form 

Imx(r,cj) = 

= ^[ln(r!,/r!„)]^-'^/'^edy„Knn(17,/17^)] ,(42) 
70W 

with 6dyn(y) = Xn{y)X^{y). After resolving the 
implicit frequency dependence brought about by the 
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moment in the definition of il^j, the leading low- 
frequency behavior of the dynamic susceptibility at 
criticality is Imx [ln(17//7ow)]'^-''/'^/(7oa'). In 
the disordered Griffiths phase, we obtain Im^ ~ 
(7ou;)'^/^- [ln(fi//7ow)]i+'^/^ The dif- 
ferences from the random transverse-field Ising model 
results^^'^^ have two reasons (i) the additional frequency 
dependence hidden in the effective damping constant 7 
and (ii) the difference between a Lorentzian spectrum in 
our case and the (5-function spectrum in the Ising case. 
The local dynamic susceptibility can be found along the 
same lines, yielding the scaling form 

Imx'°'(r,w) = 

= ^[Hnj/n^)]-''^'^'e%'^[r^Hn{ni/n^)] ,(43) 

with e^°;„(y) = Xniy). At criticality, this 

leads to Imx'°'= ^ [ln(17//7ow)]-''/'^/(7oa;), and in 
the disordered Griffiths phase, we get Imx'°'^ ^ 

(7ow)^/"-V''''+'^'^(i-"^)'*/^[ln(17//7otj)]''/^ 

V. GENERALIZATIONS 
A. Non-Ohmic Dissipation 

In this section we briefly discuss how our results 
change, if we replace the Ohmic damping term {zq ~ 2) 
in the starting action (1), (2) with a non-Ohmic term 
{zo ^ 2). We are interested in the range zq = 1 to 
oo; Zq = 1 corresponds to undamped (dissipationless) 
dynamics, 1 < zq < 2 is the so-called supcr-Ohmic case 
(damping qualitatively weaker than Ohmic damping) and 
for Zo > 2, the damping is sub-Ohmic (qualitatively 
stronger than Ohmic). 

Let us first consider sub-Ohmic damping, zo > 2. In 
this case, the crucial observation is that a single clus- 
ter with sufficiently large damping constant can undergo 
a freezing or localization transition independent of the 
bulk system. In Sec. Ill A we showed that this transi- 
tion occurs when the damping constant 7 becomes larger 
than 7c = A'-^"^'^^/''°zo/{n{zo - 2)). Within the strong- 
disorder RG, the damping constant 7 renormalizes addi- 
tivcly. Thus, even for very small bare dissipation, suffi- 
ciently large and strongly damped clusters will be formed 
under the RG (as long as ^ — *■ 00 with — > 0). Once 
they are formed, their quantum dynamics freezes. Con- 
sequently, for Zq > 2 the global quantum phase transition 
is destroyed by smearing.® 

In the supcr-Ohmic case, 1 < zq < 2, the behavior is 
less exotic. To study this case, we repeat the derivation 
of the strong disorder RG recursion relations described in 
Subsection IIIB for 1 < zq < 2. As was already pointed 
out, the multiplicative form of the recursion (14) for the 
interactions J follows directly from the structure of sec- 
ond order perturbation theory and does not depend on 
Zq. In contrast, the recursion for the gaps e does depend 



on the value of zg. Repeating the exact solution of the 
two-site cluster for the super-Ohmic case, we find 

i~^ = a[er + er] (44) 

instead of the multiplicative form (22). Here x = {2 — 
zo)/zo and a is a constant. This form also follows from 
the fact that the damping constants add, 72 =72+73, 
together with the power-law dependence e ~ ^2^0/(^0-2) 
^-1/2: of the single-cluster gap on the damping constant 
derived in (7). For undamped dynamics, zq = 1, cq. 
(44) reduces to the dirty boson result l/e2 = l/e2 + l/es 
derived by Altman et al.^^ These authors also solved the 
resulting flow equations for zq = 1 and found Kosterlitz- 
Thouless-like flows. 

While a full solution of the RG flow equations in the 
generic case 1 < zq < 2 remains a task for the future, the 
qualitative critical behavior can be inferred from the re- 
cursion relation (44). As a result of the additive form of 
(44), the local gaps e are much more weakly renormalizcd 
than the interactions which are governed by the multi- 
plicative recursion (22). Near criticality, the distribution 
of the interactions J thus becomes highly singular while 
that of the gaps e remains narrower. We therefore expect 
the critical point not to be of infinite-randomness type 
but conventional with power-law scaling t ~ , although 
the dynamical exponent z can become arbitrarily large 
as Zq — *■ 2~. Similar behavior was found at a percolation 
quantum phase transition. 

B. Generic > 1 

So far, all of our explicit calculations have been for the 
large- limit of the 0{N) order parameter field theory. 
In this subsection we show that the results do not change 
qualitatively for all A'^ > 1, i.e., all continuous symmetry 
cases. In order to do so, we reanalyze the recursion rela- 
tions (14) and (22) for generic A^ (the relations (15) and 
(21) trivially carry over for all A^). For definiteness, we 
focus on the case of Ohmic dissipation. 

As discussed above, the multiplicative form of the re- 
cursion (14) for the interactions relics on the structure 
of second order perturbation theory only, it is thus valid 
for all A^ including the discrete Ising case. In contrast, 
the form of the recursion (22), which describes how the 
local gap e (i.e., distance from criticality) changes if two 
clusters are combined, potentially does depend on A^. To 
understand this dependence, we first look at the related 
problem of the dependence of e on the size (moment) of 
the cluster. 

By invoking the quantum-to-classical mapping it was 
recently shown^^ that the gap depends exponentially on 
the size, e ~ e~^^ (with c being a constant) for all con- 
tinuous symmetry cases A^ > 1. This follows from the 
fact that classical one-dimensional continuous-symmetry 
0{N) models with 1/r^ interaction arc known to be ex- 
actly at their lower critical dimcnsion^^'^®'^^ implying an 
exponential dependence of the correlation length on the 
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coupling strength. Alternatively, one can explicitly es- 
timate the strength of the transverse fluctuations in a 
putative ordered phase and notice the logarithmic diver- 
gence of /p duj/{^iij). The exponential size (moment) de- 
pendence of the gap e requires a multiplicative structure 
of the recursion relation (22) for the merging of two clus- 
ters because their moments simply add, see Eq. (15). We 
thus conclude that this multiplicative structure is valid 
for all continuous symmetry cases, A'^ > 1. 

Consequently, for sufficiently broad disorder distribu- 
tions, the complete set of recursion relations (14), (15), 
(21) and (22) is vahd for all > 1, and with it the 
resulting infinite-randomness scenario of Sec. IIIC. Pos- 
sible A^-dependent prefactors modify nonuniversal quan- 
tities only. An analogous conclusion was drawn in the 
undamped case, zq = 1, in Refs. 53,58. 

The universal behavior of all continuous symmetry 
cases has to be contrasted with the case of Ising sym- 
metry, A^ = 1. In the Ising case, the gap does not de- 
pend exponentially on the cluster size. Instead, for suf- 
ficiently large Ohmic dissipation, the cluster dynamics 
freezes, i.e., it undergoes the localization transition of the 
dissipative two-state system. The resulting behavior of 
an Ising system with Ohmic dissipation is thus very sim- 
ilar to that of a continuous-symmetry system with sub- 
Ohmic dissipation (as discussed in the last subsection): 
Sufficiently large clusters freeze independently from the 
rest of the system which leads to a destruction of the 
global quantum phase transition by smearing. This be- 
havior was predicted in Ref. 6 and recently confirmed by 
an analytical strong-disorder RG*" as well as numerical 
simulations. 



VI. CONCLUSIONS 

In summary, we have studied quantum phase transi- 
tions in systems with continuous-symmetry 0{N) order 
parameters under the influence of both quenched disor- 
der and dissipative dynamics. To this end, we have ap- 
plied a strong-disorder RG to the LGW order parameter 
field theory of the transition. For Ohmic dissipation, we 
have found an exotic infinite-randomness critical point 
in the same universality class as the random transverse- 
field Ising chain. In the sub-Ohmic case, the quantum 
phase transition is destroyed by smearing, while super- 
Ohmic damping (including the undamped case) leads 
to conventional behavior. These results must be con- 
trasted with the case of Ising symmetry for which an 
infinite-randomness critical point occurs in the absence 
of damping'*'^ while Ohmic dissipation causes a smeared 
quantum phase transition. ^'^ 

All these different behaviors and their relations can be 
understood with the help of a general classification*'^^ of 
phase transitions in the presence of weak disorder. This 
classification is based on the effective dimensionality of 
the defects or, equivalently, the rare regions: If finite-size 
regions are exactly at the lower critical dimension of the 



problem, the critical point is of infinite-randomness type 
(accompanied by power-law quantum Griffiths singular- 
ities). Here this applies to continuous- symmetry order 
parameters with Ohmic dissipation as well as dissipa- 
tionless Ising order parameters. If the rare regions are 
below the lower critical dimension, the behavior is con- 
ventional (continuous symmetry order parameters with 
super-Ohmic dissipation); and if they are above the lower 
critical dimension, individual regions order (freeze) inde- 
pendently, leading to a smeared transition (continuous 
symmetry order parameters with sub-Ohmic damping or 
Ising systems with at least Ohmic damping). 

It is worth noting that Del Maestro et al.^^ very re- 
cently studied the large- A^ action (3) in one dimension by 
numerically solving the saddle-point equations. All their 
results are in beautiful agreement with our predictions, 
i.e., they confirmed that the quantum critical point is of 
infinite-randomness type and in the universality class of 
the random transverse-field Ising model. 

We now turn to potential experimental realizations 
of our theory. One application is the Hertz-Millis 
theory^*'^^ of the (incommensurate) itinerant antifer- 
romagnetic quantum phase transition. In this theory, 
the LGW free energy (1) is derived from a microscopic 
Hamiltonian of interacting electrons by integrating out 
the fermionic degrees of freedom in favor of the order 
parameter field Lp. While this procedure involves inte- 
grating out soft (gapless) particle-hole excitations and 
is thus potentially dangerous, '^^ the resulting order pa- 
rameter field theory of the antiferromagnetic transition 
appears to be internally consistent and free of additional 
singularities, at least in three dimensions. However, the 
applicability of the theory to realistic systems is still a 
controversial question, in particular for the much-studied 
heavy fermion compounds where several experimental re- 
sults are in pronounced disagreement with the theoretical 
predictions. ^^'^^ Different scenarios to explain the dis- 
crepancies are discussed in the literature (see Ref. 62 for 
a recent review) , and one much-discussed reason are dis- 
order effects. 

Our theory provides explicit results on how the inter- 
play of dissipation and disorder in the vicinity of the itin- 
erant antiferromagnetic quantum phase transition can 
yield activated dynamics, quantum Griffiths phenomena, 
and non-Fermi liquid behavior. We expect this to make 
an experimental verification or falsification of the disor- 
der scenario much easier. Note that a generic metallic 
system will have extra complications not contained in 
the LGW free energy (1). Specifically, attention must be 
paid to the long-range Ruderman-Kittel-Kasuya-Yosida 
(RKKY) part of the interaction between the magnetic 
fluctuations. It can produce an extra subohmic dissipa- 
tion of locally ordered clusters^^ which leads to freez- 
ing into a "cluster glass" phase at a low non-universal 
temperature Tqq determined by the strength of the 
RKKY interactions. This phase replaces part of the 
quantum Griffiths regions. It's properties and the zero 
and finite-temperature transitions to the surrounding 
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by our infinite-randomness critical point. Possible phase 
diagram scenarios are sketched in Fig. 2. 

Another potential application that has attracted con- 
siderable attention recently is the superconductor-metal 
quantum phase transition occurring as a function of wire 
thickness in ultrathin nanowires."^" The clean version of 
this transition was studied by means of a one-dimensional 
LGW theory (1) with a complex order parameter (equiv- 
alent to 2) and Ohmic dissipation. ■^^'■^^ However, 
there is experimental evidence for the pair breaking 
in this system being caused by magnetic impurities at 
the surface of the nanowire. This inevitably introduces 
quenched disorder due to the random positions of the 
magnetic impurities. Our theory thus describes the ther- 
modynamics of this quantum phase transition. With 
proper modifications, it should also apply to arrays of 
resistively shunted Josephson junctions. 

So far, we have focused on the thermodynamics close 
to the quantum phase transition. Transport properties 
can also be calculated within the strong-disorder RG by 
following the approach of Refs. 51,52. Calculations along 
these lines are underway; their results will be reported 
elsewhere. 



FIG. 2: (Color online) Schematic phase diagrams of a disor- 
dered itinerant quantum antiferromagnet close to the quan- 
tum phase transition showing two possible scenarios for the 
appearance of a cluster glass phase, which is denoted by CG. 
The other labels are as in Fig. 1. 

phases are not fully explored, yet (the transitions may be 
of fluctuation-driven first order at low temperatures^'^). 
The behavior of observables in the broad quantum criti- 
cal region above the cluster glass phase will be controlled 
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